Amplitude and Gradient Scattering in Waveguides with Corrugated Surfaces 
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We study chaotic properties of eigenstates for periodic quasi-lD waveguides with "regular" and 
"random" surfaces. Main attention is paid to the role of the so-called "gradient scattering" which is 

■ due to large gradients in the scattering walls. We demonstrate numerically and explain theoretically 

that the gradient scattering can be quite strong even if the amplitude of scattering profiles is very 

, small in comparison with the width of waveguides. 
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During last decade much attention has been paid to the theory of quasi- ID disordered solids with the so-called bulk 
scattering. By this term one describes the situation where the whole volume of a scattering region contains scatters 
whose density determines the mean free path A for a propagation of electrons. According to the theory apart from 
A, transport properties of finite samples are described by two other characteristic lengths: size L of a sample and 
C$ ■ localization length . The latter is deterimied by the degree of the decrease of the amplitude of eigenstates along 
infinite samples with the same scattering characteristics. The core of the modern theory of the transport for such 
quasi-lD systems is the so-called single-parameter scaling. It was shown that when the mean free path is much less 
than both L and Zoo, all statistical characteristics of the transport are fully described by the only scaling parameter 
which is the ratio of the localization length to the size of a sample (see, e.g. [1] and references therein). 

Another kind of quasi-lD systems that has attracted much attention in past few years, is the many-mode waveguide 
l— ~" '. with rough surfaces. In this case the scattering is entirely related to statistical characteristics of scattering walls, 
' therefore, one can speak about surface scattering. For some time it was believed that the surface scattering can 
^ . be analytically described by modified methods thoroughly developed for bulk scattering. However, recent numerical 
studies of such systems [2,3] have revealed a principal difference between surface and bulk scattering (see discussion 
and references in [4]). Specifically, it was found that the transport through quasi-lD waveguides with rough surfaces 
essentially depends on many characteristic lengths, not on one length as in the case of the bulk scattering. This 
fact is due to a non-isotropic character of scattering in the channel space. In particular, the transmission coefficient 
£f) ' smoothly decreases with an increase of the angle of incoming waves, since characterictic lengths for backscattering 
. are different for different channels [3,5]. 

The latter subject of the surface scattering has a direct link to the problem of guantum chaos. The point is that 
the waveguides with rough walls can be treated from the viewpoint of classical and quantum mechanics that describe 
a particle moving inside billiards and having multiple reflections from the walls. One of the problems of quantum 
. chaos is the quantum- classical correspondence for the situation when, in the classical limit, global properties of the 
motion of a particle are strongly chaotic. More specifically, it is of great interest to find what is the fingerprint of 
classical chaos in quantum eigenstates of closed/periodic billiards, as well as the relation of statistical properties of 
J-J ' the transport through open billiards to the underlying classical chaos. 

^ . In this paper we investigate the properties of quantum eigenstates of billiards with regular and rough walls, with the 
application to the wave scattering through quasi-lD waveguides with surface scattering. To be specific, we consider 
quasi-lD waveguides that are periodic in the x— direction with period 2k. The upper and low walls are given by 
the functions, fi(x) = fi(x + 27r) = d + ai£ i i(x) and fi(x) = f%{x + 2tt) = 02^2(2^) where d is the average width 
of the waveguides and 121,2 <C d stands for the amplitude of the scattering walls. Our interest is in the structure of 
eigenstates of the corresponding Hamiltonian with zero boundary conditions on the two walls. 

For our purpose it is convenient to pass to the variables u = x, v — fo(J)-Ji{x) * n wn i° n the new Hamiltonian 
H = H° + V describes a particle moving inside a waveguide with flat boundaries in new u, v coordinates (see details 
in [6]). Here H° = <^r(-P« + ^v) an d the effective potential V(u,v,P u ,P v ) depends on the functions fi and fa, 
with P u , P v as the canonical momentums. The solution of the Schrodinger equation for H can be written in the form 
iPe{u,v) — exp(ix) ip x (u, v). Since statistical properties of eigenstates of H do not depend on specific value of the 
Bloch index \ inside the first Brillouin band, all numerical data were obtained for a specific value of \- By expanding 
il>x( u i v ) m the basis of Ho, one can find the matrix representaion of H in the unperturbed basis specified by the two 
indexes n and m [7]. 
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The Hamiltonian matrix elements are given by 
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with /{ 2 = df li2 (u)/du. 

One should note that the matrix elements in the new variables u, v depends both on /i i2 and on their derivatives 
f[ 2 - This very fact demonstrates a highly non-trivial role of scattering walls since it is a problem to separate the 
influence of the amplitude scattering from the gradient scattering. As was shown in Ref . [8] , the perturbation approach 
is quite complicated and should be performed in a special way. 

Our goal is two-fold. First, we would like to compare the case of a "regular" upper profile = cos(w) that has 

been studied in detail in [6,7] and [9], with the "random" one £i(u) — J2n T cos(Nu) where the amplitudes Ajq 
are chosen at random and Nt — 100. In both these cases we assume ^(u) = and we refer to this case as to the 
" symmetric case" . In this way we explore the role of roughness of the surface scattering. 

Second, we wish to analyse the role of the gradient scattering that is due to the derivatives /{ 2 - Fot this we 
consider the "asymmetric case" with fi=d + a£(u) and fa = a^(u) with a = a\ = a 2 , = As one can 

see from the expressions for the matrix elements, in this case the scattering is only due to the "gradient terms" since 
fi — fa = d = const. Specifically, we have J^, n = J^, n — J®, n — 0, and the rest of the integrals are reduced to 

Jn'n - c J q ^ n '- n ^du ; J 2 , n - e 2 e u e i{n '- n)u du ; J 3 n , n = ^ J e^'-^du . (5) 
where = d£(x)/du, and e = a/d. 

One natural representation of the Hamiltonian matrix Hi^i{\) =< I \ H \ V > x is the "channel representation" 
for which one fixes the values of n starting from the lowest one, n — —N maXl and running over all values of m. In 
numerical simulations we have to make a cutoff for the values of m and n in the Hamiltonian matrix. Our data refer to 
the ranges, 1 < m < M max and \n\ < N max with N max — 32 , M max — 62, for which the total size of the Hamiltonian 
matrix is L = (2N max + l)M max = 4030. 

^,From the structure of the Hamiltonian matrices shown in Fig.l one can make some important conclusions. First, by 
passing from "regular" wavegudes with Nt = 1 to "random" ones with Nt = 100 in both symmetric and asymmetric 
cases, the Hamiltonian matrices tend to be fully filled by off-diagonal elements. However, the matrices correponding 
to the random waveguides keep the block structure, thus indicating some regularity in spite of a completely random 
character of the surface scattering. This fact was shown to result in a kind of " non-ergodicity" for the eigenstates of 
the total Hamiltonian H for any high energy, see details in Ref. [10]. It is clear that standard theoretical approaches 
based on completely random matrices can not adequately discribe the structure of eigenstates. As was shown in Ref. 
[5], for open waveguides with rough surfaces and large number of channels there are many characteristic lengths, in 
contrast to the bulk scattering where the so-called single-parameter scaling holds. 

First, Figure 1 clearly demonstrates that the eigenstates of H are expected to be much more extended in the 
unperturbed basis for rough profiles, than for one-cosine profiles. Second, it is instructive to make a comparison 
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between the symmetric and asymmetric cases (compare (a) with (b) and (c) with (d) in Fig.l). The data show 
that the matrices for the assymetric cases have smaller elements than those for the symmetric ones. This is in 
correspondence with the fact that large part of terms in the expressions for the matrix elements of H vanishes due 
to the absence of the amplitude scattering. Therefore, one can expect that the eigenstates are less random for the 
asymmetric case. 




FIG. 1. Central part of the Hamiltonian matrix. The 62 x 62 blocks corresponding to n, n' = [—3, 3] are shown. The larger the 
amplitudes of the matrix elements the more black the corresponding regions are. Four cases refer to: (a) one-cosine symmetric 
waveguide, (b) one-cosine asymmetric waveguide, (c) symmetric waveguide with a rough wall, Nt = 100, (d) asymmetric 
waveguide with rough walls, Nt = 100. 



In order to analyze the structure of eigenstates of H in detail, we have diagonalizcd the Hamiltonian matrices 
shown in Fig.l and constructed the "state matrices" IC 1 "! 2 . Here Cf are the amplitudes of the eigenstates in the 
basis representaion given by the index I. Namely, the index / refers to unperturbed basis states that correspond to 
the unperturbed Hamiltonian H° . The index a refers to a specific exact eigenstate. All eigenstes are reordered in 
increasing energy, with a = 1 the ground state. Therefore, to understand how strongly localized/extended are the 
exact eigenstates in the unperturbed basis, one should fix the value of a and explore the dependence of |C"| 2 on I. 

Let us now compare the symmetric cases with the asymmetric ones, see Fig. 2. The most important conclusion 
that can be deduced from the data is that the eigenstates are, in general, more extended in the asymmetric cases. 
Specifically, for the same small values of a there arc more components C" with large values of I, in comparison 
with the symmetric case, compare (c) with (d). At a first glance, this looks strange since for asymmetric profiles 
the Hamiltonian matrices are obviously less "random" than for the symmetric ones, as is mentioned above. Close 
inspection of the data in Fig. 2 shows that there is an additional effect that is also important in connection with the 
structure of eigenstates. Namely, the eigenstates for asymmetrical cases (b) and (d) turn out to be more sparse in 
comparison with the cases (a) and (c). This is manifested by a large number of "holes" along each line in (b) and 
(d) for fixed values of a, in the comparison with the cases (a) and (c). Thus, the eigenstates for waveguides with 
only gradient scattering (asymmetric case) are more extended in the unperturbed basis, and, at the same time, more 
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sparse than for the wavegudies with both gradient and amplitude scattering (symmetric case). 

This phenomenon is important in view of the scattering properties through waveguides of finite size with the profiles 
such as considered here. As is known, chaotic structure of eigenstates of closed (or periodic) waveguides/billiards is 
directly related to the scattering properties of open systems with the same profiles. For example, the degree of local- 
ization of eigenstates of closed systems determines the degree of localization of scattering states, and correspondingly, 
the value of transmission through open systems. 
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FIG. 2. Lower part of the state matrix | | 2 . The data for all cases (a-d) refer to the same profiles for which the structure 
of Hamiltonian matrices is shown in Fig.l. Here darker regions correspond to larger values of |C™| 2 . 



In conclusion, we have studied the structure of eigenstates of quasi- ID periodic waveguides with regular and random 
walls. Main attention was paid to the role of gradient scattering in comparison with the amplitude scattering. It was 
shown that for the case when the scattering walls have large number of harmonics, the gradient scattering itself is very 
strong. This is demonstrated by the data obtained for the waveguides with asymmetric walls, for which the amplitude 
scattering is absent. It was also revealed that the role of the gradient scattering is highly non-trivial. Specifically, the 
gradient scattering turns out to be relatively strong in the absence of the amplitude scattering. 

We are grateful to N. Makarov for fruitful discussions. This work was supported by the CONACyT (Mexico) Grant 
No. 34668-E, and IIG3G02, VIEP, BUAP. 
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